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We study a one-dimensional SU(N) Hubbard model with an attractive on-site interaction and 
N > 2 at half-filling on the bipartite lattice using density-matrix renormalization-group method 
and a perturbation theory. We find that the ground state of the SU(N) Hubbard model is a charge 
density wave state with two-fold degeneracy. All the excitations are found to be gapful, resulting in 
an insulating ground state, on contrary to that in the SU(2) case. Moreover, the charge gap is equal 
to the Cooperon gap, which behaves as —2Nt 2 /(N — 1)U in the strong coupling regime. However, 
the spin gap A s and the quasiparticle gap Ai as well open exponentially in the weak coupling 
region, while in the strong coupling region, they linearly depend on U such that A s ~ —U(N — 1) 
and Ai U(N -l)/2. 



I. INTRODUCTION 

Correlation effects in electronic systems have been of 
long-term interest in the condensed matter physics. In 
recent years, important progress has been made experi- 
mentally in ultra-cold atomic systems where strong corre- 
lation leads to some novel physical phenomena. In par- 
ticular, interacting fermionic atoms can be trapped in 
an optical lattice^. More interestingly, the interaction 
in the ultra-cold atomic systems is tunable through the 
Feshbach resonance, which allows for a full exploration 
of various fundamental properties of strongly correlated 
models. Moreover, the nuclear spin of the atoms can 
be larger than the electronic spin. One could expect 
richer physics induced by the degree of freedom of higher 
spins^££. Therefore, one would like naturally to gen- 
eralize the Hubbard model with two spin components to 
the one with iV-components. 

In this paper, we investigate low energy properties of a 
one-dimensional half-filled SU(N) Hubbard model 7 - with 
an attractive on-site interaction and N > 2. This is 
a generalization of our previous work in which we have 
focused on the SU(4) case^. The Hamiltonian of the one- 
dimensional SU(N) Hubbard model is represented by 

l u L 

n = -t^2 'YjA r c i+la + h.c.) + — n^i*' > 

i=l <r i=l a ^ a ' 

(1) 

where t > is the hopping matrix set as the energy unit, 
U < the coupling constant of the attractive on-site in- 
teraction, L the number of lattice sites, a and a the spin 
indices, which take the values (N - l)/2, (N - l)/2 - 1, 
• • •, 1 - (N - l)/2, -{N - l)/2, respectively. c\ a and 
Cjo- denote the creation and annihilation operators, re- 
spectively, for a particle with spin a at the site i and 
JT-ia = cl a Ci a indicates the operator of the particle num- 
ber. The Hamiltonian (1) has U(1)®SU(N) symmetry^. 
The generator for the U(l) symmetry is J\f = ^2 ia 



and the U(l) symmetry implies that the particle number 
TV is conserved. The generator for the SU(N) symmetry 
is S A = J2iaa' Hrr^~ A '^ia' ' wnere T A are the generators 
of the SU(N) group in its fundamental representation. 
These symmetries are useful for simplifying numerical 
calculations and classifying excitations for the present 
system. 

Although this model is exactly solvable for N = 3», 
it seems that physical features obtained for N = 2 are 
not immediately applicable to more general cases with 
N >2H. Recently, this model with N > 2 and U > has 
been studied by using several analytic as well as numer- 
ical approaches. At half-filling, a renormalization-group 
analysis 7 shows that both the charge and spin coupling 
constants are renormalized to a large value, resulting in 
gaps in both sectors. An analytic perturbative renor- 
malization group treatment in the fermionic representa- 
tion alternatively gives rise to a gapful spectrum for all 
N > 2 and U > 0^2. Moreover, by employing bosoniza- 
tion method and quantum Monte Carlo simulations, As- 
saraf et al found for a 1/N filling that the spin excitation 
is gapless for any positive U, whereas the charge excita- 
tion is gapful only for U > U c where U c ^ 0. However, 
Buchta et al obtained gapless spin as well as gapful charge 
excitations from accurate density-matrix renormalization 
group(DMRG) calculations with N = 3,4 and 5 for any 
U >0. 

For the attractive interaction, on the other hand, it 
is known that the one-dimensional attractive half-filled 
SU(2) Hubbard model is described by a Luther-Emery 
liquid model, in which the charge excitation is gapless, 
whereas the spin excitation is gapful. By the hidden 
SU(2) transformation, the SU(2) Hubbard model with 
U can be mapped to the one with —U. However, for 
the SU(N) Hubbard model with N > 2 such a mapping 
does not exist so that one cannot obtain any insight into 
the low-energy properties for the attractive case through 
the mapping from the repulsive case. In this paper, we 



2 



will show that the SU(N) Hubbard model at half filling 
with an attractive interaction belongs to a different uni- 
versality class from the SU(2) one and all the excitations 
are gapful. We expect that our findings are not only 
of fundamental interest, but also useful for experimen- 
talists, since the attractive SU(N) Hubbard model may 
be possibly realized by experiments in ultra-cold atomic 
systems. 

The rest of this paper is organized as follows. In Sec. 
II, we analyze low-energy properties of the Hamiltonian 
(TI]) in the strong coupling limit (JV — 1)|£^| 3> t first 
and then in the weak coupling limit (N — l)\U\ <C t by 
the perturbation treatments. Particularly we discuss the 
dependence of charge gap, Cooperon gap, spin gap and 
quasiparticle gap on U in both regimes. In Sec. Ill, we 
present our numerical results which arc obtained from 
the DMRG calculations and compare them with analytic 
behavior in both weak and strong coupling regime. A 
summary is finally given in Sec. IV. 



II. PERTURBATION CALCULATIONS 

In this section, we study the low-energy properties of 
the Hamiltonian by a perturbation theory. For this 
purpose, we rewrite it as 



TL = Tit + Ti v 



(2) 



where Tit = ~tJ2ia(Ci<jC i + la ' + h.c.) is the hopping term 
and 7i u — [U /2)J2ia^ia' ^io-^ia' ^ ne on_ site interaction. 
We start with the on-site interaction part TL U . Since the 
on-site interaction is attractive, N particles with differ- 
ent a tend to stay on one site and form a SU(N) sin- 
glet. The energy of the SU(N) singlet is UN(N - l)/2. 
In the half-filling case where N/2 particles per site, the 
ground states of TL U are highly degenerate, involving half 
of the lattice sites occupied by the SU(N) singlets and 
the other half being empty. In the strong coupling re- 
gion, i.e. \U\(N — 1) 3> t, Ti u is taken as the zeroth order 
Hamiltonian, while Tit, being the order of Nt, is regarded 
as a perturbation. Up to the second order, we obtain the 
first order effective Hamiltonian 

H$ f = PH t P = 0, 

where P is a projection operator which restricts the effec- 
tive Hamiltonian in the subspace spanned by the ground 
states of H u - The second order effective Hamiltonian 
reads 



H e/j = Pn t ——-(\-p)n t p 

2t 2 



1 



-—p y^hiP 

1)U ^ 

% 



(N-l) 
2t 2 



where hj = Y]^ hi a the number operator at the site i. 
The degeneracy of the ground state of Ti u is thus lifted by 
Tit- At half filling, one can easily obtain the ground state 
energy correction per site which is given as Nt 2 / (N — 1)U 
by the first term of Eq. $3§. Since the second term in- 
troduces an effective repulsion interaction between the 
particles at nearest neighbor(NN) sites, it results in a 
charge density wave (CDW) state in which the SU(N) 
singlet and the empty site occur alternatively to exhibit 
a long range order resulting from translational symme- 
try breaking. It turns out that the ground state is two- 
fold degenerate. The configuration of the ground state is 
schematically shown in Fig. QJa). 
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FIG. 1: Schematic illustrations of the configurations, and the 
empty circle here represents one SU(N) singlet formed by N 
particles with different a. (a) For the CDW ground state. 
Shifting all the circles by one site, one can obtain the other 
one of two-fold degenerate ground states, (b) For the charge 
excitation: the SU(N) singlet at site i + 2 shifts to the site 
i + (c) For the Cooperon excitation: one SU(N) singlet is 
added to the site i+1. 

It is well-known that both the charge and Cooperon 
excitations for the attractive SU(2) Hubbard model are 
gapless. However for the SU(4) case, it has been shown 8 
that the charge and Cooperon excitations are gapful and 
equal to each other at half-filling. In the following, one 
can see that this conclusion is also valid for other attrac- 
tive SU(N) Hubbard models with N > 2. In particular, 
the charge gap and Cooperon gap can be easily derived 
from the effective Hamiltonian ©. The charge gap A c 
is defined as the lowest excitation in the SU(N) singlet 
subspace as follows 



A c = Ei(L, NL/2, 0) - E (L, NL/2, 0), 



(4) 



where Eq(L, M, S) is the ground state energy in the spin- 
s' channel with L sites and M particles, and E n (L, M, S) 
the n-th excitation energy. The gap for Cooperon exci- 
tations for the SU(N) Hubbard model is given from the 
energy difference between states by adding N particles 
or iV holes to the system, which is defined as 



, NL 



N,0) 



, NL 
ML,— 



N,0) 



NL 

-E (L,—,Q). 



(5) 
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The charge gap A c is given by shifting one of SU(N) 
singlets in the ground state configuration to its nearest 
neighbor site as shown in Fig. QJb). Then one has 



2Nt 2 



(N-1)U 

Similarly, one obtains the SU(N) Cooperon gap 

2Nt 2 



(N — 1)U' 



(0) 



(7) 



which is shown in Fig. \V[c). One can see that A c = 
for all N > 2, which was previous shown for the N = 4 
case&. In the large- N limit, 



Ar, = A 



N 



2Nt 2 2t 2 



We note that the difference for the low-lying excitations 
between N = 2 and N > 2 results from the effective in- 
teraction between the singlets at the NN sites as involved 
in effective Hamiltonian For the N = 2 case, one has 
the effective Hamiltonian 14 



H 



(2) 

e//,s«(2) 



9/ 2 



(9) 



Compared to the effective Hamiltonian ([3]), one has ad- 
ditionally a pair hopping term, which involves the same 
amplitude as the NN repulsion term and eventually de- 
stroys the CDW long range order for N = 2. On the 
other hand, for N > 2 cases, although a similar hopping 
term emerges at the 7V-th order perturbation, it has a 
smaller amplitude than the NN repulsion term so that 
one can has a stable CDW ground state. 

Now we turn to study the spin and quasiparticle ex- 
citations. The spin gap is defined in correspondence to 
the lowest excitation with different spin quantum num- 
ber from the ground state. The quasiparticle gap is de- 
fined as a energy change by adding one particle or hole 
to the system. Since the ground state is CDW, in or- 
der to obtain these two gaps, we resort to the following 
Hartree-Fock(HF) approximation. 



n ir ,n-i ~ n. 



'{ fl ia ') + {n i(T )n ia ' - (rW)(n iCT >), (10) 



where {hi a ) = no + (— l) l Sn, and Sn is the order param- 
eter, (■ • ■) is the average over the ground state. The HF 
Hamiltonian then reads 



n HF = -*^(4c i+lCT + /i. c .) + 

+ {n%o)n ia > - {h i<T )(n i(7 <}). 



U 



(ii) 



At half-filling, one has no = \ and ^ Zcr fii a = Eq. 
(fTTj) can be further simplified as 

H HF = -tJ2(4a^+la + h.C.) 

+ U(N - l)5nJ2{~l) l n ia + const. (12) 



Introducing ai a = cW, bi a 
Fourier transformation &i a = 

T J2k bkoe tkl , we obtain 



c-2i+io and taking a 



H HF = -tJ2((l+e~ ik )*l*h*+h.c.) (13) 

ha 

+ U(N - l)Sn ^{a\ a a k(T - b\J> ka ) + const. 



Diagonalizing this HF Hamiltonian, we can obtain two 
bands for each spin species with the following dispersion 



± 



Af+4< 2 cos 2 ^, 



(14) 



where A i = —(N — l)USn is the quasiparticle gap. Under 
the HF approximation, one has that A c = A s = 2Ax- By 
solving the self-consistent equation for Sn = {aj a ai a — 
I! S| 5 ) one can determine the order parameter Sn. In 
the weak coupling limit t <C |f/|(JV — 1), we can solve it 
approximately and obtain 



Ai ~ 27rie Lr <«- 1 ) . 



(15) 



In the strong coupling limit, one can put Sn ~ 0.5, and 
then obtain the quasiparticle gap 



Ai ~ —U(N - l)/2, 



and the spin gap 



A., 



-U(N-1). 



(16) 



(17) 



In the HF approximation, since A c = 2Ai, one obtains 
A c ~ ~U(N — 1) which is inconsistent with the results 
Eq. ([5]) obtained from the strong coupling perturbation 
theory. This is because in the presence of the strong at- 
tractive on-site interaction between particles, the single 
particle picture is not valid for the charge excitations, 
in which N particles forming a SU(N) singlet shift as a 
whole to the NN site, as illustrated in Fig. [U(b). How- 
ever, we would expect that the HF results are still qual- 
itatively correct for weak coupling limit. 



III. NUMERICAL CALCULATIONS 

In this section, we present our numerical results for 
N = 3, 4, 5 and 6 cases and compare them with those an- 
alytic predictions from the perturbation theories and the 
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HF approximation. Our numerical results are obtained 
from a large scale DMRG computation ! 15 ' 16 ! 17 It is well 
known that DMRG method is the most powerful numeri- 
cal tool for accurate exploration on low-energy properties 
of one dimensional systems at zero temperature. How- 
ever, it is nontrivial to apply this method to the SU(N) 
Hubbard model. On one hand, there is a large number of 
degrees of freedom per site. For instance, the degree of 
freedom per site is 32 when N = 5. On the other hand, 
although a open boundary condition (OBC) in DMRG 
calculations may lead to more accurate results than a 
periodic boundary condition (PBC), while the emergence 
of edge states under the OBC makes the calculation of 
various excitations practically more difficult. 

For these reasons, we have made some additional con- 
siderations to the standard DMRG algorithm as follows. 
When N = 3, i.e. SU(3) case, one has 8 degrees of free- 
dom per site. In this case, we use the PBC with keeping 
about 2000 ~ 3500 optimal states in the DMRG pro- 
cedure and two sites are added in order to enlarge the 
chain length at each DMRG step. The maximal trunca- 
tion error is of the order 10 -5 . Similarly, the charge gap 
for N ~ 4 is calculated with at most 1600 states kept. 
The other excitations are calculated with the OBC. At 
each step we add one lattice site to the chain, break- 
ing the lattice site into two pseudo sites. To obtain ac- 
curate results, we preselect some specific chain lengths, 
and perform sweeping at these preselected lengths, the 
final results are got by extrapolating the results at the 
preselected lengths to the thermodynamic limit. To avoid 
cumbersome edge excitations when the OBC is employed, 
the preselected lengths are odd instead of even. Corre- 
spondingly the definitions of the excitations are changed. 
For example, the Cooperon gap for N — 4 is redefined as 
A N = E (L,2L + 6,0)-E (L,2L + 2,0)-W, where the 
particle-hole symmetry is taken into account and odd L 
are used. 



A. Energy correction and degeneracy of the 
ground states 

The analysis based on the strong coupling perturba- 
tion theory in the above section shows that the hopping 
term, lifting the degeneracy of H u , results in an energy 
correction per site Nt 2 /(N — 1)U. To verify this predic- 
tion, we calculate the energy corrections for N — 3,4 and 
5 using DMRG method. Both numerical and analytic re- 
sults are shown in Fig. [5] One can see that the energy 
correction is negative and decreases monotonically with 
respect to increasing of U. In the strong coupling region, 
the results given by the perturbation theory agree well 
with the DMRG data. For N = 3, the deviation between 
analytic and numerical results is within the numerical 
accuracy for U < —2. With increasing of N, the devia- 
tion becomes smaller and smaller. This implies that the 
perturbation theory gives rise to better results for larger 
N for the strong coupling regime. In the weak coupling 




FIG. 2: Ground state energy correction per site obtained from 
both the DMRG calculation and the strong coupling pertur- 

(2) 

bation theory 7~i e ff for N=3, N=4 and N=5 are shown as a 
function of U. 

region, the analytic results from the effective Hamilto- 
nian ([3]) severely deviates from the DMRG results. This 
is reasonable since in the weak coupling region, TL t is 
no longer a small perturbation. Nevertheless, when N 
increases, the valid region of the strong coupling pertur- 
bation theory increases simultaneously. 
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FIG. 3: Two lowest excitation gaps in the SU(N) singlet sub- 
space for N = 3 case with two different values of U. In the 
thermodynamic limit, the first excitation energy becomes de- 
generate with the ground state, while the energy difference 
between the second excitation and the ground state remains 
finite, indicating a gapful excitation. 

We have examined the symmetry properties of the 
ground state. In the thermodynamic limit, the ground 
state is two- fold degeneracy. For a finite L, the ground 
state is a SU(N) singlet, and there is another SU(N) sin- 
glet state just above the ground state. When the chain 
length L is increased, the energy difference between these 
two states decreases and eventually vanishes in the large- 
L limit. This is consistent with the prediction from the 
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effective Hamiltonian ([3]). Moreover, for a given chain 
length L, we find that the energy difference between two 
states for the Hamiltonian |T]) decrease rapidly with in- 
creasing \U\, since the effective Hamiltonian Eq. ^ be- 
comes more accurate for larger |Z7|. In Fig. [3] and [4] we 
demonstrate this feature with two different values of U 
under OBC for N = 3 and 4, respectively. Note that the 
energy difference for the second excitation, which is fi- 
nite in the large- L limit, corresponds to edge excitations 
rather than bulk excitations. 
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FIG. 4: Fwo lowest excitation gaps in the SU(N) singlet sub- 
space for N = 4 case with two different values of U. In the 
thermodynamic limit, the first excitation energy becomes de- 
generate with the ground state, while the energy difference 
between the second excitation and the ground states remains 
finite, indicating a gapful excitation. 



B. Charge gaps and Cooperon gaps 

The prediction that the charge gap A c is equal to the 
Cooperon gap Ajv is confirmed numerically for N=4 8 . 
However, the strong coupling perturbation theory gives 
rise to A c = A^r for general N, it is still interesting to ex- 
amine numerically whether this equation is valid for both 
even and odd N and all regimes of the coupling constant 
U— . Fig. [5] demonstrates the charge and Cooperon gaps 
for N — 3 as comparison to N = 4 as a function of U. 
One can see that both charge and Cooperon gaps for both 
N = 3 and N = 4 are finite for all U < 0. When \U\ is 
increased, these gaps increase first in the weak coupling 
region. After they reach their maxima, they decrease fol- 
lowing the asymptotic behavior — nfjrfwf resulting from 
the strong coupling perturbation theory. We would men- 
tion that even beyond the valid range of the perturbation 
theory, our numerical data still show that the charge gaps 
are equal to the corresponding Cooperon gaps for both 
odd N(= 3) and even N(— 4). Therefore, one may con- 
clude that A c = Ajv for all N and U < 0. 

Moreover we calculate the Cooperon gaps for N = 5 
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FIG. 5: The Cooperon and charge gaps are shown for N = 3 
and 4. The asymptotic behavior from the effective Hamilto- 
nian© (curves) are also shown for comparison. 



and 6 in order to explore the lavge-N behavior of the 
charge as well as Cooperon gaps. The results of the 
Cooperon gap for N = 5 and 6 are shown in Fig. [6] 
Comparing Fig. [5] and El one can see that the larger N 
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FIG. 6: The Cooperon gaps and the corresponding asymp- 
totic behaviors as a function of U for N=5 and 6 are shown. 

is, the closer Cooperon gaps to that derived from the ef- 
fective Hamiltonian ([3]). For N = 3 as shown in Fig. 
the deviation of the asymptotic behavior from the DMRG 
data is visible even at U = —10, while for N = 4 case, the 
deviation is within the numerical accuracy up to U ~ —5. 
For TV = 5 and 6, the good agreement can be seen up to 
U ~ —3 and U ~ —2, respectively. Furthermore, as N 
increases, the height of the peak of the Cooperon gaps in- 
creases and the position of the peak shifts toward U = 0. 
These features can be easily understood since the crite- 
ria for the strong coupling perturbation t <C |C/|(-ZV — 1) 
depends on |U|(iV - 1) rather than U only. The DMRG 
results here also verify that the charge and Cooperon gap 
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TABLE I: Decomposition of some irreducible representations 
of SU(3) to SO(3). 
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behaves asymptotically as 
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in the large- TV limit. 
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C. Quasiparticle gap and spin gap 

In the following, we discuss the quasiparticle and spin 
gaps. In order to improve numerical efficiency, we need 
to carry out DMRG calculations in a subspace which 
can be obtained by decomposing the irreducible repre- 
sentations of SU(N) into the irreducible representation 
of SO(3) 18 ' 19 . Numerically, one can determine the ir- 
reducible representations for short chains by varying z- 
component. Table U shows the decomposition of some 
irreducible representations of SU(3) to SO (3). 



TABLE II: Decomposition of some irreducible representations 
of SU(4) to SO(3). 



SU(4) 


SO(3) 


V 


[I 4 ] 





1 


[2 2 ] 


4029290 


20 


[2 1 ! 2 ] 


3©2©1 


15 


[4 1 ] 


0©2©3©4©6 


35 


[3 1 ! 1 ] 


Iffilffi2ffi3ffi3©4©5 


45 


[I 1 ] 


3 
2 


4 


[2 2 ! 1 ] 


2 W 2 fcB 2 fcB 2 


20 


[3 1 ! 2 ] 




36 



For N — that the spin excitation belongs to the 
representation [2 1 ! 1 ] and it is 8-fold degenerate, while 
the quasiparticle excitation belongs to the representation 
[l 1 ] and it is 3-fold degenerate. Table HT1 shows the de- 
composition of some irreducible representations of SU(4) 
to SO (3). In this case, the spin excitation belongs to 
the representation [2 1 ! 2 ] and its degeneracy is 15-fold, 
whereas the quasiparticle excitation belongs to the rep- 
resentation [l 1 ] and its degeneracy 4-fold. Fig. [7] shows 
the quasiparticle and spin gaps for N = 3 and 4 as well. 
In weak coupling region, one can see that both the quasi- 
particle and spin gaps open exponentially, as predicted 
by the HF approximation, although we cannot determine 
precise critical behavior for the opening of the gaps from 
our numerical data due to limited numerical accuracy. 
In the strong coupling region, both the quasiparticle gap 



FIG. 7: The quasiparticle gap and spin gap for N=3, 4 are 
shown in the figure. For N=3, the quasiparticle gap(O) and 
spin gap(D) are shown by the black symbols. For N=4, the 
quasiparticle gap(A) and spin gap(x) are shown by red sym- 
bols. 



and spin gap depend linearly on U as seen in the in- 
set. This is also consistent with the HF results. More- 
over, one also see that the ratio Ai/|[/| approaches 1 for 
N = 3 asymptotically, while it approaches 1.5 for N = 4. 
These asymptotic behaviors qualitatively coincide with 
Eq. (HU) and HH) and the relation A s = 2Ai derived 
from the HF approximation. For finite U , however, the 
relation is just approximately valid since there remains 
residual interaction between quasiparticles. 



IV. CONCLUSIONS 

In this paper, we have investigated low-energy proper- 
ties of the SU(N)(N>2) Hubbard model with attractive 
on-site interaction at half-filling. By using the perturba- 
tion theory and DMRG method, we show that the ground 
state is a CDW state with two-fold degeneracy. The 
CDW long range order reflects the broken translational 
symmetry. On one hand, the strong coupling perturba- 
tion theory predicts that both the charge excitations and 
Cooperon excitations are gapful and equal to each other 
with the asymptotic behavior —2Nt 2 /(N — 1)U. Com- 
bining this with numerical results for N — 3,4, 5 and 6, 
we conclude that the charge gap is equal to the Cooperon 
gap for all U < and N > 2. In the large N limit, both 
the charge and Cooperon gaps behaves as —2t 2 /U. Con- 
sidering the CDW ground state, we obtain qualitatively 
correct behavior for the spin and quasiparticle gaps from 
HF approximation, which are confirmed by our numerical 
data. In the weak coupling region, they open exponen- 
tially, whereas in the strong coupling region, they depend 
linearly on U with the spin gap A s ~ — (N — 1)U and 
quasiparticle gap Ai ~ — (N — l)U/2. 

Our findings indicate that at half filling, the 
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SU(N)(N>2) Hubbard model belongs to a different uni- 
versality class from the SU(2) case. This can be easily 
understood by observing the difference between the ef- 
fective Hamiltonians for the SU(N)(N>2) and the SU(2) 
cases. In particular, the effective Hamiltonian which for 
the SU(2) case involves an effective hopping term has 
the same order as the effective repulsion interaction for 



particles between nearest-neighbor sites. The hopping 
term destroys the CDW order. For N > 2, the repulsive 
interaction dominates over the effective hopping term. 
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